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We report the measurement of the transverse momentum dependence of inclusive J/il) polarization 
in p + p collisions at ^ = 200 GeV performed by the PHENIX Experiment at RHIC. The J/ip 
polarization is studied in the helicity, Gottfried- Jackson, and Collins-Soper frames for pt < 5 GeV/c 
and \y\ < 0.35. The polarization in the helicity and Gottfried- Jackson frames is consistent with 
zero for all transverse momenta, with a slight (1.8 sigma) trend towards longitudinal polarization 
for transverse momenta above 2 GeV/c. No conclusion is allowed due to the limited acceptance 
in the Collins-Soper frame and the uncertainties of the current data. The results are compared 
to observations for other collision systems and center of mass energies and to different quarkonia 
production models. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

Quarkonia production in high-energy hadronic colli- 
sions is an essential tool for investigating QCD. The QQ 
pair is produced in a hard scattering involving gluons, 
which is followed by a hadronization process that forms 
the bound state. These formation and hadronization 
steps are the subject of many studies. Initial tests of 
quarkonia production models using J/ip cross sections 
measurements are still inconclusive fl!|, suggesting that 
other observables would be useful to challenge the dif- 
ferent production models. For example, a key piece of 
information to help pin down the mechanism of heavy 
quarkonia (cc and 66) production and the bound state 
formation is the angular distribution of its decay leptons. 

The angular distribution of spin-i lepton decay from 
quarkonium (spin 1) is derived from density matrix ele- 
ments of the production amplitude and parity conserva- 
tion rules [2-4]. The angular distribution integrated over 
the azimuthal angle is given by 



da 



dcosO" 



= A(l + Acos2 r) 



(1) 



where A is a normalization factor and 6* is the angle 
between the momentum vector of one lepton in the po- 
larization quarkonium rest frame and the longitudinal 
direction (z coordinate) of a selected polarization vector 
(frame). The polarization parameter A is related to the 
diagonal elements of the density matrix of the produc- 
tion amplitude and contains both the longitudinal (ctl) 
and transverse (ctt) components of the quarkonium cross 
section 



A = 



aT — 2a-i, 
ctt + 2aL 



(2) 



The quarkonium polarization is longitudinal (transverse) 
in a given frame if A is negative (positive) . 

The most common polarization frame used in analy- 
ses performed at collider experiments is where z is the 
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FIG. 1: Definition of the polarization frames: helicity(HX), 
Gottfried- Jackson(GJ) and Collins-Soper(CS) frames. 



quarkonium momentum. Polarization measured in this 
manner is referred to as being in the helicity frame 
(HX) 0. In fixed target experiments the most fre- 
quently used polarization frame has z as one of the col- 
liding hadrons momentum in the quarkonium rest frame, 
namely, the Gottfried- Jackson frame (GJ) 0. Another 
polarization frame, used primarily for the studies of 
Drell-Yan production, is the Collins-Soper frame (CS) Q 
that defines z as the bisector between the directions of 
the first colliding parton and of the opposite of the sec- 
ond colliding parton in the dilepton rest frame. A dia- 
gram representing the three polarization frames is shown 
in Fig. [H The amplitude and the sign of A depend on 
the frame used in the measurement. The natural polar- 
ization axis for the production process can be defined as 
that where the lepton decay azimuthal angle distribution 
is symmetric and A is maximum In such a frame, the 
density matrix of the production amplitude is diagonal. 

Several quarkonium production models have been pro- 
posed to describe the perturbative terms which are rele- 
vant for QQ production, while other models include non- 
perturbative terms related to the formation of the bound 
state. The various models predict different polarizations 
and are described below. 

In the Color Evaporation Model (GEM) quarko- 
nia production is assumed to be a fixed fraction of the 
pQCD cross section for invariant masses between twice 
the mass of the heavy quark (c or h) and twice the mass 
of the open heavy quark meson [D or B). This model has 
reasonable agreement with most of the measured quarko- 
nia cross sections but no predictive power for the polar- 
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ization (3]. Nevertheless, according to [sl], multiple soft 
giuon exchanges destroy the polarization of the heavy 
quark pair. 

The earliest Color Singlet Model (CSM) was a calcula- 
tion of the leading order gg ^ S-wave charmonium + g 
process where the relative momentum of the QQ pair 
with respect to the quark mass mn is neglected and 
the pair is produced on-shell |9l-[l]|. The QQ bind- 
ing is calculated from potential model wave functions. 
J/V' yield measurements reported by CDF [l2[ and 
PHENIX [l] are largely underestimated by this model. 
The J/ip polarization predicted by LO CSM is transverse 
in the HX frame fl3| . Subsequent calculations also in- 
cluded NLO terms nO-Iil, NNLO terms [13, d and an 
s-channel cut contribution that allows off-shell cc quarks 
to end up in the bound state J^j- These calculations 
show large changes in the yield and polarization relative 
to the earlier calculations. The new calculations of the 
J/ip yield is closer to what is observed in PHENIX and 
CDF for pt < 10 GeV/c. The J/ip polarization is pre- 
dominantly longitudinal in the HX frame according to 
these new calculations. 

Non-relativistic QCD (NRQCD) effective theory ^ 
makes use of short distance (rng) and non-relativistic 
(mQiy'^^ terms, where is the typical quark velocity in 
the quarkonium rest frame. A typical v for charm (bot- 
tom) is 0.3c (0.1c). The S-wave charmonium is described 
as a series of intermediate color singlet^^^^ or color octet'-^-' 
state contributions 



(8) 

J 9 



(3) 



+ 0{.')\'s['^gg)+0{.')\'S^'^g 



where the spectroscopic notation is used. The 

non-perturbative operators O (v) are parametrized us- 
ing experimental results. Since the singlet state has a 
small contribution to the yield, this model is also referred 
as the Color Octet Model (COM). Using constraints 
from the CDF cross section J/ip data, reasonable agree- 
ment is obtained with PHENIX yield results assuming 
J/ ip production is dominated by gluon fusion in the ^S'q^^ 
and ^Pq^^ intermediate states for pT < 5 GeV/c [T3l.[2l|. 
Calculations performed in estimated X(^Sq^^) = 
and X{'^Pq^^) = —0.05 indicating a very small longitudi- 
nal polarization from direct J/ip in this px range. Nu- 
merical estimations [l^, [13] and subsequent NLO correc- 
tions [23 supports that the polarization for pr ^ Mj/,^, 
where production from gluon fragmentation is supposed 
to be important, is predominantly transverse in the HX 
frame. 

The J/ip polarization in hadronic collisions was stud- 
ied in fixed target experiments at y^s < 39 GeV 26-33]. 
These experiments predominantly covered \xf\ > and 
Pt < 5 GeV/c. In @ it was noted that J/ip polarization 
measured in the CS frame by HERA-B [31| and by 
E866/NuSea [S^ smoothly changes from longitudinal to 
transverse with the total momentum. The polarization 



observed in CDF at ^/s = 1.96 TeV in midrapidity for 
Pt > 5 GeV/c showed a small longitudinal polarization 
in the HX frame [34| . This result contradicts the first LO 
CSM and COM expectations. 

Complimentary J/ip polarization measurements in 
p+p collisions at ^/s — 200 GeV can help elucidate the 
production mechanism. Moreover, it is expected that 
the polarization of J/ ip is modified in the presence of nu- 
clear matter effects in d+Au collisions and hot and dense 
matter in Au-|-Au collisions j35| . Thus future measure- 
ments of J/ ip polarization in Au and Au-I- Au at RHIC 
demands a good reference from p+p collisions. 

This paper reports the transverse momentum depen- 
dence of the J/ip polarization for \y\ < 0.35 in the HX, 
GJ and CS reference frames. The study was performed 
in the dielectron decay channel for pT < 5 GeV/c . The 
experimental apparatus used to measure electron decays 
from J/ip mesons is detailed in section [Hi The proce- 
dure followed to obtain the cos 6* distributions, the cor- 
responding polarization parameters and their uncertain- 
ties are explained in section iHll The results, comparison 
with measurements at other facilities and interpretation 
in the context of current theoretical models are presented 
in section HVl 



II. 



EXPERIMENTAL APPARATUS AND 
J/iP IDENTIFICATION 



This analysis was performed with data collected in 
p+p collisions at = 200 GeV during the 2006 RHIC 
Run with the PHENIX central arm detectors [36]. The 
geometrical coverage for single electrons corresponds to 
pseudorapidity ]77] <G.35. Each one of the roughly back- 
to-back two arms covers A(p = tt/2. 

Data were recorded using a minimum-bias trigger that 
required at least one hit in each of the two beam-beam 
counters (BBC) located at 3.0 < \r]\ < 3.9 and scanning 
approximately 50% of the p+p cross section. A dedicated 
trigger (EMCal RICH Trigger - ERT) was also used to 
select events with at least one electron candidate. The 
ERT required a minimum energy in any 2x2 group of 
the Electromagnetic Calorimeter (EMCal) towers ^ and 
associated hits in the Ring Imaging Cerenkov detector 
(RICH) in coincidence with the minimum-bias trigger 
condition. The EMCal energy threshold was set to 0.4 
GeV and 0.6 GeV for two different periods during the 
data taking run. 

Collisions within ±30 cm of the center of the detector 
along the beam direction were used in this analysis. After 
data quality selection, the number of collisions sampled 
was 143 billion BBC triggers, corresponding to an inte- 
grated luminosity of J £ = (6.2 ± 0.6) pb~^. 



^ Corresponding to A'q X Acp = 0.02 X 0.02 rad 
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FIG. 2: Invariant mass of dielectrons in the J/V' mass range. 
Dashed Unes represents the mass range used in the polariza- 
tion analysis. 



Electron candidates were selected from tracks recon- 
structed in the Drift Chamber (DCH) and in the Pad 
Chamber (PC) with momentum larger than 0.5 GeV/c. 
Electron identification was achieved by requiring the 
tracks to be associated with at least one fired phototube 
within a ring radius 3.4 cm < Rrmg < 8.4 cm centered 
on the projected track position in the RICH. In addi- 
tion, the presence of a matching energy cluster in the 
EMCal was required within four sigma in both the po- 
sition and expected energy/momentum ratio. Since the 
hadronic background in the J/ip mass region is small in 
p+p collisions, only loose electron identification criteria 
were used. 

Dielectron pairs from J/ijj decays were counted in 
the invariant mass range € [2.9,3.2] GeV/c^. The com- 
binatorial background was estimated using like-sign 
(e+e"*" and e~e~) pairs. Since we evaluated the ERT ef- 
ficiency using J/ ip simulation, we required that the ERT 
segment was fired by one of the J/tp decayed electrons. 
Hence, only pairs with at least one electron matching 
geometrically the position of an actual ERT trigger in 
the event were accepted. The dielectron mass distri- 
bution in the J/^p mass region is shown in Figure [2] 
The signal/ (combinatorial background) ratio was 28. Af- 
ter combinatorial background subtraction, we counted 
2442 ± 51 e+e" pairs withp^ < 5 GeV/c in the selected 
J/ip mass range. These counts include a residual con- 
tinuum background, which consists mainly of correlated 
open heavy quark decays to electrons. This background 
was found to be less than 10%. 



III. ANALYSIS PROCEDURE 

The angular dependence of the detector response for 
electrons from J/tp decays was estimated using a full 
GEANT3 [131 based detector simulation. Dead chan- 
nels or malfunctioning regions were removed from the 
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FIG. 3: Distribution of single electrons versus the Z coordi- 
nate of the track in the drift chamber for real (open boxes) 
and simulated (closed boxes) data for different sectors (0 is 
the bottom and 3 is the top ones) in the east (left) and west 
(right) detector arms. The box-height for each point corre- 
sponds to its statistical uncertainty. 



detector simulation and from the real data analysis. The 
experimental acceptance was checked by simulating sin- 
gle electrons with collision z vertex and px distributions 
weighted to reproduce observed distributions of electron 
candidates from real data. The simulated detector accep- 
tance for these single electrons was compared to that for 
real data for different azimuthal sectors (see Fig. [3]). Re- 
maining differences in the acceptance between simulated 
and real data were attributed to conversions 7 — > e~^e~ 
in the detector support structure at large z which were 
not included in the simulation. These differences were 
accounted for in the systematic uncertainty listed in Ta- 
ble H 

The detector response to electrons in the simulation 
was tuned to match the data. A clean sample of elec- 
trons in the data was obtained by selecting electrons from 
Dalitz decays and photon conversions in the beam pipe 
which were identified by their very low invariant mass 
[ssj . Figure |4] shows the comparison between the momen- 
tum dependent electron identification efficiency (ee/u) 
for single electrons from fully reconstructed Dalitz and 
photon conversion decays in minimum-bias data and sim- 
ulation. Good agreement above 0.5 GeV/c was achieved 
within the statistical uncertainties. 

The pt dependence of the ERT efficiency was esti- 
mated for each one of the 8 EJVICal sectors by taking 
the fraction of minimum-bias single electron candidates 
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TABLE I: Systematic uncertainties in the pr dependent polarization measurement in the hehcity and Gottfried- Jackson (in 
parentheses) frames. 



Description 


[0, 1] GeV/c 


[1,2] GeV/c 


[2, 5] GeV/c 


[0, 5] GeV/c 


Acceptance 


0.006 (0.036) 


0.006 (0.012) 


0.006 (0.008) 


0.006 (0.024) 


Polarization bias in acceptance 


0.022 (0.047) 


0.0011 (0.005) 


0.008 (0.031) 


0.012 (0.032) 


Continuum fraction 


+0.033 (+omi\ 
-0.021 \-0.0U) 


+0.023 /'+0.032\ 
-0.027 (,-0.062j 


+0.014 /'+0.023\ 
-0.039 \-0.070) 


+0.019 /+0.032\ 
-0.026 l-0.058j 


Input pt in simulation 


0.034 (0.062) 


0.005 (0.049) 


0.024 (0.028) 


0.034 (0.054) 


Input y, Z vertex in simulation 


0.000 (0.007) 


0.000 (0.007) 


0.000 (0.007) 


0.000 (0.007) 


Run-by-run fluctuations 


0.019 (0.123) 


0.016 (0.035) 


0.016 (0.020) 


0.017 (0.050) 


ERT efficiency 


0.017 (0.110) 


0.015 (0.051) 


0.018 (0.024) 


0.015 (0.043) 


TOTAL 


+0.06 l+Q.21\ 
-0.05 \-0.W) 


+0.03 /+0.09\ 
-0.04 \-0.W) 


0.04 (!?,;««) 


0.05 (in?) 



W 

1 

0.8 
0.6 
0.4 
0.2 




Dalitz + y ^ e"^e" decays 




O real data 
■ simulation 



1 2 3 4 5 

electron momentum [GeV/c] 

FIG. 4: Single electron identification efhciency estimated us- 
ing full reconstructed Dalitz decays from real data (open cir- 
cles) and simulation (full squares). Dotted line represents the 
minimum pr for the electron used to reconstruct J/^ decays. 



X! 

^10 



-a 



-a 



10- 



10- E 



p+p^J/\|/ + X lykO.35 \/^ = 200GeV 



da 



dydp^ [l+[v^lhf\^ 



8 9 10 
[GeV/c] 



FIG. 5: Fit to J/^ yield times dielectron branching ratio [B) 
after detector acceptance and efficiency corrections for the 
real data with A = 28.7 ± 1.0 nb/GeV/c, h = 3.41 ± 0.21 
GeV/c, and n = 4.6 ± 0.4. 



(Fig. [5]) , and a Gaussian function was fit to the ra- 
pidity dependence of the J/ if) yield reported in [l[ 
and to the collision Z vertex distribution. 



that fired the ERT. These efficiencies were used in the 
ERT simulation. Changes in the trigger thresholds and 
channel masks in the ERT during the run period were 
used in the simulation in order to reproduce realistic run 
conditions. 

The tuned detector simulation was used to reproduce 
the measurement of J/t/i dielectron pairs and to match 
their momentum, rapidity and vertex distributions. The 
kinematics of the simulated J / ip were estimated in four 
steps: 

1. Unpolarized J/ip e~^e~ pairs were generated with 
uniform distributions in rapidity (|y| < 0.5), 
Pt {pt < 7 ^^), azimuthal angle (— tt < < tt), 
and collision vertex along the beam axis Z 

{\Zvertex\ < 40 Cm). 

2. The J/ip Pt distribution obtained after applying 
the efliciency and acceptance corrections agrees 
with the previous result A Kaplan function 
__da_ _ Apr ^g^g ^Yie pj, distribution 



The fitted px , rapidity and collision vertex 
functions were then used to re-weight the simu- 
lated J/ip events. The top half of each plot in 
Figs. [6] and [7] shows the cos 61* distributions in 
the HX, GJ, and CS frames of e+e~ pairs in the 
J/ijj mass range obtained in J/ijj simulation and 
real data ^ after combinatorial background sub- 
traction. The simulated and real data distribu- 
tions are functions of the detector acceptance and 
efficiency and the original dN ^+ ^- / dcos 9* in the 
J /ip mass range. The bottom panels show the ratio 
between the real data and simulated A = distri- 
butions, corresponding to the acceptance corrected 
cosO* distributions. 



dydpi 



\l+{PT/by 



The cos 8* resolution estimated in the simulation was 0.08 
in the HX, 0.025 in the GJ and 0.007 in the CS frames. 
These resolutions are much smaller than the bin width of the 
cos 6* distributions used in the polarization analysis. 
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o 
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0<p^<lGeV/c GJ 
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l<p^<2 GeV/c 



0.5 1 -1 
cos 9* 



-0.5 



0.5 1 

cos 9* 



2<p^<5 GeV/c 



-0.5 
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cos 9* 




0.5 1 1 

cos 9* 



0.5 1 
cos 9* 



FIG. 6: Each of the six plots shows (upper half) cos 5* distributions of positrons decayed from (solid points) measured and 
(horizontal bars) simulated J/^p mesons and (lower half) acceptance corrected distributions obtained from the ratio between 
real and simulated J/tp distributions. Fits to Eqn. 1 are represented as solid lines. Dashed lines correspond to one standard 
variation of the parameters in the fit. The top three plots are for the helicity (HX) frame and the bottom three are for the 
Gottfried- Jackson (GJ) frame, where the smaller cos 0* range causes larger uncertainties of the fits. 



Equation ([T]) was fitted to these acceptance cor- 
rected 008 6** distributions with no constraints on 
the parameters. Solid lines are the most likely fits 
and dashed lines represent 68% confidence level in- 



terval. In the CS frame, the fit returned a po- 
larization which was out of the physical limits 
(A S [—1,1]). This was a result of the small ac- 
ceptance for the cos 6* distribution in the PHENIX 
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FIG. 7: (top) Same as in As in Fig. [6] but now for the CoUins- 
Soper frame. Dotted lines on bottom panel correspond to 
Eq. d) for A = ±1. 



central arms for this frame, leading to a large statis- 
tical uncertainty on its polarization measurements. 
Thus, the CS frame is no longer considered in this 
article. 

4. Any asymmetry in the electron decay distribution, 
i.e. A 7^ 0, can change the detector acceptance. 
Hence, the fourth and final step of the simula- 
tion was to apply a weight in cos6'* to the sim- 
ulated J/ip by using the A obtained in the third 
step. When using this realistic angular distribu- 
tion for the pt dependent acceptance, and the cor- 
responding uncertainties, we obtained a variation 
in the yield up to ±8% for px < 5 GeV/c that cor- 
responds to changes in polarization results no larger 
than 0.02 in the HX frame and 0.05 in the GJ frame. 
These variations were accounted for in the system- 
atic uncertainties. 

We also estimated the contribution to the J/ tp polari- 
zation from the continuum background by measuring A 
in the dielectron mass range [1.7,2.3] GeV/c^. The ac- 
ceptance and efficiency corrections were performed using 
simulated DD — > e+e~ decays, the dominant source of 
e+e" pairs in [1.7, 2.3] GeV/c^, according to the analysis 
in [39} . The polarization in this mass range is consistent 
with zero, with values between ±0.3 in the HX and ±0.9 
in the GJ frame. The 10% continuum contribution can 
change the measured polarization in the J/Tp mass range 
by at most 1qq2 HX frame and 14 in GJ frame and 
was included in the systematic uncertainties. 

The A measurement is also sensitive to differences be- 
tween acceptance in simulated and in real data, run-by- 



TABLE II: J/t/j polarization results in the helicity and 
Gottfried- Jackson frames. Transverse momentum is in 
GeV/c. Uncertainties correspond to statistical and system- 
atics respectively. 



Pt 



(PT ) 



A 



A 



0- 1 

1- 2 

2- 5 
0-5 



0.64 
1.47 
2.85 
1.78 



PI -I r-l-0.12 +0.06 
'^•-'-^-0.18 -0.05 



_n in+o-09 +0-03 

"•^'-'-0.13 -0.04 



-o.i9i;]:i^ ±0.04 



-0.10lo:og ±0.05 



|-| f.1 +0.39 +0.21 
'-'■'-'^-0.52 -0.19 



-0.20 



+0.30 +0.09 
.0.32 -0.10 



-0.35 



+0.18 +0.06 
-0.22 -0.09 



n -ie+0.18 +0.09 
U.ro_g j2 -0.11 



run condition variations, uncertainties in rapidity, Z ver- 
tex, and transverse momentum shape inputs to the simu- 
lation, as well as the ERT efficiency px shape. These un- 
certainties were introduced as variations in the efficiency 
and weighting parameters for different detector sectors in 
the simulation. Resulting variations in A were accounted 
for as systematic uncertainties and are listed in Table 
m The systematic uncertainties are correlated between 
different px ranges. The total systematic uncertainty 
is taken to be the quadratic sum of these components, 
assuming they are uncorrelated. Additional checks in- 
cluded the variation of the minimum momentum require- 
ment of the single electrons and the rejection of tracks 
going to the edges of the detector. These variations re- 
turned only statistical fluctuations in the polarization re- 
sults. 



IV. RESULTS 

Figure [H] shows the transverse momentum dependence 
of the J/tp polarization in the HX and JG frames. The 
uncertainties of the fit are larger in the GJ frame given 
the smaller cos 9* range compared to HX frame. The nu- 
merical values are listed in Table [Hi For the HX frame 
also shown are currently available theoretical models: 
COM [H and the s-channel cut GSM ^ calculated us- 
ing the same polarization frame. There are no theoretical 
predictions for the GJ frame. 

The measurements presented here are for inclusive 
J/V' ■ Feed-down from Xc and ip' may also contribute 
to the observed polarization and are not separated out. 
The world average result for the feed-down contribution 
to the J/ip yield is 33 ± 5% [40| . The polarization of the 
indirect J/ip should be smeared during the decay pro- 
cess. If the J/ip from feed-down sources are unpolarized, 
the direct J/ip may have a larger A in magnitude than 
that reported here. 

The J/V' polarization is consistent with zero for all 
transverse momenta but exhibits a 1.8 sigma longitudinal 
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FIG. 8: J/ijj polarization parameter (Aj/^) versus transverse 
momentum (pr). Boxes are correlated systematic uncertain- 
ties, (upper) Helicity frame data is compared with COM[i3] 
and s-channel CSM [3| calculated in the same polarization 
frame, but there is no prediction for CEM. (lower) There are 
no theoretical predictions for the Gottfried- Jackson frame. 



FIG. 9: xf dependence of J/'0 polarization for pt < 
5 GeV/c measured by PHENIX, HERA-B [SH], E771 ^ 
and E866/NuSea (CS frame) [11]. 



polarization at pT > 2 GeV/c in the HX and GJ frame 
when the quadratic sum of the statistical and system- 
atic uncertainties are considered. In the HX frame, the 
PT dependent A follows the s-channel cut CSM expecta- 
tions for prompt J/^p fiTl. Finally, the COM prediction 
[Ta] . using the NRQCD matrix elements fitted to CDF 
data, is also consistent with our data over the px range 
covered by the calculation. 

Figure|9]shows that the polarization for < 5 GeV/c 
follows what is observed in fixed target experiments for 
a more extended xp range in the HX and GJ frames. 
Statistical and systematic uncertainties are quadratically 
summed for this comparison. Note that the E866/NuSea 
result was measured in the CS frame. 



In principle, intermediate singlet and octet color states 
may be absorbed differently in the nuclear matter present 
for fixed target p+A measurements, possibly changing 
the final J/tp polarization. However, comparisons are 
limited by the uncertainties in the present data. The 
observed agreement between the p + p results reported 
here and fixed target p + A measurements are not yet 
able to determine the magnitude of nuclear matter ef- 
fects on J/ip polarization. Direct comparison between 
future high statistics pt and rapidity dependence of the 
J/ip polarization in p+p and o?-|-Au collisions will provide 
a better picture for these effects. 
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V. CONCLUSIONS 

We have presented the first J/ij} polarization measure- 
ment at RHIC for two different polarization frames. The 
observed pt -dependent J/?/) polarization parameter in 
the HX frame is consistent with the s-channel cut CSM, 
COM and no polarization within current uncertainties. 
The integrated momentum polarization observed in both 
the HX and GJ frames are in good agreement with the 
results obtained at fixed target experiments collected in 
lower energy p+A collision in the same xf region. Up- 
coming higher luminosity p+p data will allow more ac- 
curate measurements over the full decay angular distri- 
butions and over extended pt and rapidity ranges. 
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